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HIGHLIGHTS 


•  Three  mouse  strains  were  subjected  to  a  stress  model  eliciting  PTSD-like  features. 

•  Under  stress,  three  strains  demonstrated  different  phenotypic/genetic  aspects. 

•  BALB/cj  displayed  delayed  arousals  of  many  phenotypic  features. 

•  C57BL/6j  showed  signatures  of  poor  cardiac  health. 

•  Genetic  heterogeneity  and  stress  history  contributed  in  PTSD  prevalence. 
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Systematically  distinguishing  genetic  liability  from  other  contributing  factors  is  critical  for  designing  a 
preventive  strategy  for  post-traumatic  stress  disorder  (PTSD).  To  address  this  issue,  we  investigated  a 
murine  model  exposing  C57BL/6j,  DBA/2j  and  BALB/cj  mice  to  repeated  stress  via  exposure  to  conspecific 
aggressors  (Agg-E).  Naive  mice  from  each  strain  were  subjected  to  the  proximity  of  aggressor  (Agg)  mice 
for  6  h  using  a  ‘cage-within-a-cage’  paradigm,  which  was  repeated  for  5  or  1 0  days  with  intermittent  and 
unpredictable  direct  contact  with  Agg  mice.  During  the  Agg-E  stress,  DBA/2j  developed  a  different  strategy 
to  evade  Agg  mice,  which  potentially  contributed  to  its  phenotypic  resilience  to  Agg-E  stress.  Although 
Agg  mice  inflicted  C57BL/6j  and  BALB/cj  with  equivalent  numbers  of  strikes,  BALB/cj  displayed  a  distinct 
behavioral  phenotype  with  delayed  exhibition  of  a  number  of  PTSD-like  features.  By  contrast,  C57BL/6j 
mice  displayed  unique  vulnerability  to  Agg-E  stress  induced  myocardopathy,  possibly  attributable  to  their 
particular  susceptibility  to  hypoxia.  A  group  of  genes  {Bdnf,  Ngf,  Zwint,  Cckbr,  Slc6a4,  Fkbp5)  linked  to  PTSD 
and  synaptic  plasticity  were  significantly  altered  in  C57BL/6j  and  BALB/cj  Agg-E  mice.  Contributions  of 
Agg-E  stress  history  and  genotypic  heterogeneity  emerged  as  the  key  mediators  of  PTSD-like  features. 
Linking  genetic  components  to  specific  phenotypic  and  pathological  features  could  have  potential  clinical 
implications. 

©  2015  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Trauma  exposure  is  a  necessary,  but  not  sufficient  condition  for 
developing  PTSD,  which  may  explain  why  only  a  sub-population  of 
trauma  survivors  develops  PTSD.  The  prevalence  of  PTSD  appears 
to  be  mediated  by  genetic  predisposition  putatively  vulnerable  to 
the  environmental  factors  [1,2].  The  gene-environment  interaction 
model  explains  the  higher  PTSD  prevalence  among  monozygotic 
twins  as  compared  to  their  dizygotic  counterparts  exposed  to 
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equivalent  traumatic  burdens  [3-5].  Similarly,  higher  PTSD  preva¬ 
lence  was  reported  among  the  offspring  of  holocaust  survivors 
[6]  and  of  families  with  troubled  histories  [7],  In  the  context  of 
accumulating  evidence  suggesting  that  PTSD  could  be  a  high-to- 
moderately  heritable  disease,  molecular  markers  independent  of 
environmental  effects  could  have  significant  clinical  importance 
[8-10],  However,  finding  appropriate  cohort  of  controls  [1,8]  and 
controlling  for  magnitude  of  traumatic  burden  remain  major  chal¬ 
lenges,  as  environmental  interference  can  potentially  overwhelm 
genetic  liability  [11], 

Addressing  these  challenges  requires  in  vivo  investigations 
[12],  Inter-strain  genetic  heterogeneity  and  the  histological  and 
endocrinological  variables  are  often  considered  to  be  the  media¬ 
tor  of  the  strain-to-strain  endophenotypic  variances  in  mice.  The 
inter-strain  variability  reflected  in  the  sizes  of  the  hippocampus 
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[13],  nucleus  accumbens  [14]  and  neocortical  elements  [15],  as 
well  as  the  diversities  in  the  cholinergic  [16]  and  dopaminergic 
[17]  feedbacks  could  contribute  to  social  and  cognitive  plasticity 
[18,19], 

Hence,  the  present  study  is  focused  on  characterizing  the  con¬ 
tributions  of  genetic  heterogeneity  and  environmental  challenges 
in  developing  PTSD-like  features  elicited  by  an  agonistic  social 
situation  including  both  threat  and  attack.  As  developed  by  our 
group  [20],  the  Agg-E  repeated  stress  model  attempts  to  emu¬ 
late  the  combat-like  situation,  where  soldiers  are  exposed  to  grave 
threats  repeatedly  and  uncontrollably.  Arguably  combat-related 
PTSD  symptoms  could  differ  from  those  of  civilian  PTSD  patients. 
Our  model  concurred  with  the  etiological  and  face  validators  man¬ 
dated  for  a  PTSD  rodent  model  [12].  Using  C57BL/6j  mice  as 
Agg-E  cohort,  this  model  characterized  the  genomics  [21,22]  and 
metabolomics  [23]  signatures  of  PTSD.  The  present  model  utilizes 
SJL  male  mice  trained  to  be  aggressive  toward  the  intruders;  hence 
unlike  the  typical  social  stress  models  [24,25],  the  present  model 
limited  the  direct  intruder-resident  interactions  to  not  more  than 
one  minute;  longer  direct  interaction  could  cause  serious  and  unde¬ 
sired  injury  in  Agg-E  mice.  Such  aggressive  interaction  is  designed 
to  make  the  Agg-E  mice  encounters  as  repeated  life-threatening 
events,  a  critical  feature  of  a  combat  theater.  Typical  PTSD  models, 
such  as  the  brief  footshock  model  [14,26]  and  the  predator  stress 
model  [27]  fails  to  emulate  certain  key  features  of  combat  theater 
such  as  the  repeated,  random  and  uncontrolled  occurrence  of  the 
potentially  life-threatening  events  and  the  conspecific  interactions, 
respectively. 

Nevertheless,  all  of  these  models  demonstrated  the  inter¬ 
play  of  gene  and  environment  in  developing  fear  responses  and 
other  key  PTSD-like  features.  Both  the  magnitude  and  character¬ 
istics  of  the  traumatic  load  play  a  major  role  in  determining  the 
strain-to-strain  differences  in  the  stress  adaptation  strategies.  Eor 
instance,  the  social  stress  model  elicited  acute  social  avoidance  in 
BALB/cj  in  comparison  to  C57BL/6j  [25],  The  single  foot  shock  aug¬ 
mented  social  impairment  in  DBA/2j  and  behavioral  “depression” 
in  C57BL/6j  [28]  accompanied  by  diminished  contextual  fear  con¬ 
ditioning  [14],  The  predator  stress  model  [29]  found  BALB/c  mice 
more  resilient  than  C57BL/6  and  DBA/2  strains.  With  exposure  to 
the  contextual  threat,  all  three  strains  behaved  equivalently;  but 
shifting  the  cues’  character  from  the  feline  feces  to  clay  triggered 
anxiogenic-like  responses  in  C57BL/6  and  DBA/2  mice  more  than  in 
the  BALB/c  strain.  Conversely,  a  set  of  heterogeneous  stressors  com¬ 
posed  of  open  field,  elevated  maze  and  exposure  to  novel  objects 
elicited  higher  fear  responses  in  BALB/cj  compared  to  C57BL/6j 
[30];  and  DBA/2j  performed  better  than  C57BL/6j  in  accomplishing 
passive  avoidance  tasks,  but  not  in  complex  tasks  [31], 

Motivated  from  these  studies,  we  focused  on  characterizing 
the  implications  of  the  genetic  heterogeneity  among  three  strains 
of  mice  (C57BL/6j,  DBA/2j,  and  BALB/cj)  responding  to  the  Agg- 
E  stress,  eliciting  PTSD-like  behaviors.  The  phenotypic  variability 
along  with  the  expressions  of  a  selective  genomic  panel  associated 
with  PTSD  and  synaptic  plasticity  were  evaluated. 

2.  Material  and  methods 

2.3.  Mice  obtained  for  the  study 

We  purchased  four  strains  of  male  mice,  SJL,  C57BL/6j,  DBA/2j 
and  BALB/cj,  from  Jackson  Laboratory  (Bar  Harbor,  ME),  all  at 
the  age  of  six  weeks  and  weighing  20-22  g.  Mice  were  kept  in 
a  temperature-controlled  room  (21±2“C)  on  reverse  12h/12h 
light/dark  cycle  (lights  off  at  6;00  AM)  and  all  experiments  were 
conducted  during  the  dark  phase  of  dark/light  cycle. 

SJL  albino  mice  (Agg  mice)  were  trained  to  reliably  and  quickly 
initiate  aggressive  behavior  against  the  intruder  Agg-exposed  mice 


(Agg-E  mice)  following  a  protocol  described  by  our  previous  com¬ 
munication  [20].  A  separate  room  was  dedicated  exclusively  for 
the  other  three  strains;  C57BL/6j,  DBA/2j  and  BALB/cj.  Mice  were 
single-housed  in  polycarbonate  cages  (48  cm  x  27  cm  x  20  cm)  for 
at  least  one  month  before  initiating  the  protocol. 

All  animal  experiments  were  approved  by  the  Institutional  Ani¬ 
mal  Care  and  Use  Committee  (lACUC)  at  the  Walter  Reed  Army 
Institute  of  Research,  Silver  Spring,  MD,  and  the  Medstar  Research 
Institute,  Washington  Hospital  Center,  Washington,  DC,  and  were 
performed  in  a  facility  accredited  by  the  Association  for  the  Assess¬ 
ment  and  Accreditation  of  Laboratory  Animal  Care  International 
(AAALAC). 

2.2.  Agg-E  stress  sessions 

We  described  the  detailed  protocol  in  an  earlier  communica¬ 
tion  [20],  Briefly,  the  Agg-E  stress  session  was  conducted  during 
the  dark  light  cycle.  We  caged  the  Agg-E  mouse  inside  a  small  wire 
mesh  cage  (17.5  cm  x  14  cm  x  7.5  cm),  which  was  placed  inside  the 
Agg’s  polycarbonate  home  cage  for  6  h  (9;00  AM-3 ;00  PM)  without 
any  access  to  food  or  liquid  (6h  ‘cage-within-cage’  confinement). 
The  Agg  mouse  had  nourishment  ad  libitum  and  free  access  around 
the  mesh  cage.  In  addition  to  the  indirect  exposure,  we  also  directly 
exposed  the  ‘intruder’  Agg-E  mice  to  the  ‘resident’  Agg  mice  two  to 
three  randomly  selected  times  per  day  (d ).  Each  bout  of  direct  expo¬ 
sure  lasted  for  1  min  or  until  10  strikes  were  inflicted  on  the  Agg-E 
mouse,  whichever  came  first.  In  a  separate  room,  the  control  mice 
were  kept  in  the  same  type  of  6  h  ‘cage-within-cage’  confinement 
with  fresh  bedding  for  6  h  without  any  proximity  or  exposure  to  Agg 
mice.  After  the  6  h  confinement,  the  controls  and  Agg-E  mice  were 
returned  to  their  respective  single-housing  home  cages  with  food 
and  liquid  ad  libitum  until  the  start  of  next  day’s  6  h  sessions,  which 
were  repeated  daily  for  either  5  or  10  consecutive  days,  creating 
5-day  and  10-day  stress  groups,  respectively. 

We  randomized  the  controls  and  Agg-E  mice  subjected  to  5-day 
Agg-E  sessions  into  two  groups;  one  of  the  two  groups  was  sac¬ 
rificed  1  day  after  the  5-day  Agg-E  session  and  the  other  group’s 
sacrifice  was  delayed  by  1.5  weeks  (w).  Similarly,  the  controls  and 
Agg-E  mice  subjected  to  10-day  Agg-E  sessions  were  sacrificed 
either  1  d  or  4w  post  Agg-E  sessions.  Mice  spent  their  respective 
1.5  w  or  4  w  delay  intervals  singly  housed  inside  their  home  cages 
with  food  and  water  ad  libitum. 

2.3.  Measurement  of  bodyweight,  temperature  and  urine 
markings 

The  detailed  protocol  for  recording  this  set  of  physiological 
information  was  published  in  an  earlier  report  [20].  Briefly,  body- 
weight  was  recorded  before  (Fig.  S2A)  and  after  (Fig.  lA)  each  6h 
‘cage-in-cage’  confinement  as  were  the  body  temperatures  (Fig. 
S2B-C).  The  core  body  temperatures  were  recorded  using  Elec¬ 
tronic  ID  TranspondersTM  (IPTT-300,  Bio  Medic  Data  Systems  Inc., 
DE)  subcutaneously  implanted  into  the  dorsal  cervical  region  of  the 
mice  3  d  in  advance  of  the  first  stress  session.  Territorial  distribu¬ 
tions  of  urine  marking  were  recorded  using  blotting  papers  (0.8  mm 
thickness)  placed  daily  under  the  mesh  cages  of  controls  and  Agg- 
E  mice  during  6h  ‘cage-in-cage’  confinement  (Fig.  IB).  The  papers 
were  UV-scanned  using  Molecular  Imager  Ex®  (BioRad,  Hercules, 
CA)  and  the  number  of  urination  marks  and  the  areas  contouring 
the  urine  markings  were  measured  using  Quantity  One®  software 
(BioRad,  Hercules,  CA). 

2.4.  Behavioral  assessment  during  the  direct  exposures  of  Agg-E 
to  Agg  mice 

We  visually  inspected  the  direct  interactions  between  the  Agg 
and  Agg-E  mice.  Trained  technicians  counted  and  timed  the  bouts 
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C.  Jumping  to  escape  from  aggressor  D.  Fightback  to  aggressor 


Fig.  1.  Phenotypes  of  Agg-E  mice  during  the  Agg-E  sessions.  All  records  were  taken  during  the  5-day/lO-day  Agg-ExS  (N=  16-21  mice/group);  Welch’s  t-test:  0.001, 

•p<0.05. 

In  A  and  B,  the  horizontal  dotted  line  passing  through  1.0  of  Y-axis  marks  the  baseline  defined  by  the  normalized  phenotypes  exhibited  by  the  control  mice.  MeaniSEM  of 
Agg-E  mice  divided  by  the  strain-specific  controls  (Agg-E  mice/control)  are  plotted  against  each  day  point.  (A)  Pre-stress  bodyweight.  (B)  Urine  markings. 

In  C  and  D,  behaviors  were  monitored  during  the  daily  bouts  of  direct  exposures;  the  fractions  of  the  total  Agg-E  mice  group  are  plotted  against  each  day  point.  (C)  Jumping 
to  escape  from  aggressor.  (D)  Fightback  to  aggressor. 


of  strikes,  and  monitored  three  behavioral  phenotypes  of  the  Agg-E 
mice  described  in  Table  1,  namely  jumping  (Fig.  )C),  freezing  (Fig. 
S2D)  and  fightback  (Fig.  ID). 

2.5.  Cardiac  histological  analysis 

Following  a  protocol  previously  described  [20],  heart  tissues 
collected  from  the  sacrificed  mice  were  perfused  in  ice-cold  4% 
paraformaldehyde.  Tissues  were  then  stained,  sliced,  and  mounted, 
and  a  board-certified  veterinary  pathologist,  unaware  of  the  ani¬ 
mal’s  behavioral  history  or  treatment  group,  analyzed  the  samples 
using  bright  field  optics  (Table  2). 

2.6.  Phenotypic  evaluation  during  the  partition  test 

The  partition  test  [20]  was  designed  to  present  the  contextual 
cues  reminiscent  of  the  past  agonistic  interactions  between  the  Agg 
and  Agg-E  mice.  We  performed  this  test  1  d  and  1.5  w  after  the  5- 
day  Agg-E  session  and  1  d  and  4  w  after  the  10-day  Agg-E  session. 
We  divided  the  soiled  Agg  home  cage  with  a  plastic  fenestrated 
partition  and  placed  one  Agg  and  one  Agg-E  mice  on  opposite  sides 
of  the  partition  for  5  min.  The  partition  permitted  the  passage  of 
sensory  cues  but  prevented  direct  physical  contact.  The  ethogram 
of  the  Agg-E  and  control  mice  was  recorded  and  the  videos  were 
analyzed  by  Ethovision  XT  v.7  software  (Noldus®,  Leesburg,  VA)  at 
a  rate  of  15  samples  per  second  using  dynamic  subtraction  detec¬ 
tion,  with  the  object  always  darker  than  background,  erosion  and 
dilation  filters  of  one  pixel,  and  one  sample  interval  for  averaging 
filter  [20[. 

Table  1  describes  the  behavioral  parameters  we  evaluated 
herein.  For  measurement  purposes,  we  hypothetically  divided  the 
side  of  the  cage  temporarily  holding  the  Agg-E  mouse  (or  the  con¬ 
trol)  into  three  zones  (Fig.  SI ).  The  six  inch  wide  zone  adjacent  to 


the  partition  was  defined  as  the  “partition  zone”  and  the  six  inch 
wide  region  across  the  three  side  walls  at  the  rear  was  defined  as 
the  “peripheral  zone”.  The  central  part  of  the  cage  surrounded  by 
the  partition  and  peripheral  zones  was  defined  the  “mid  zone”.  We 
defined  the  frequency  of  time  spent  in  the  partition  zone  as  the 
“partition  vigilance"  and  the  ratio  of  the  time  spent  between  the 
peripheral  and  partition  zone  as  the  “withdrawal  toward  peripheral 
region".  “Grooming"  monitored  visually  by  an  observer  is  a  com¬ 
plex  sequence  of  licking  of  paws  and  washing  of  nose,  head  and 
body.  “Retarded  locomotion"  refers  to  the  fraction  of  total  move¬ 
ment  that  happened  inside  the  partition  zone  and  “freezing"  refers 
to  an  event  of  immobility  displayed  inside  the  partition  zone  when 
all  visible  movement  of  the  body  except  the  movements  for  res¬ 
piration  became  non-existent.  “Overall  activity"  is  defined  as  the 
average  velocity  of  the  mice  recorded  during  the  5  min  test  in 
all  three  zones.  Our  previous  communication  [20]  discussed  all  of 
these  phenotypic  endpoints  in  detail.  The  behavioral  parameters 
of  the  Agg-E  mice  (Fig.  2  and  Fig.  S3)  were  normalized  by  their 
age-  and  strain-matched  controls  as  described  earlier  [20].  Present 
strategy  potentially  minimized  the  behavioral  variances  attributed 
to  the  non-Agg-E  stress-mediated  factors,  such  as  age,  food/liquid 
restriction,  movement  restriction  and  prolonged  social  isolation 
in  order  to  extract  the  unique  contributions  of  Agg-E  repeated 
stress. 

2.7.  Hemibrain  isolation  and  downstream  qPCR  analysis 

Animals  were  euthanized  by  cervical  dislocation  1  h  after  the 
partition  test.  Left  hemibrains  were  collected  from  the  mice  after 
euthanization  and  were  immediately  snap-frozen.  The  samples 
were  archived  in  a  -80  °C  freezer  for  long  term  storage  until  the  day 
of  nucleic  acid  isolation  using  Trizol  (Life  Technologies,  Carlsbad, 
CA)  following  the  vendor’s  recommendations. 
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Table  1 

Descriptions  of  the  behavioral  phenotypes  monitored  and  measured  during  the  direct  exposure  to  Agg  (Agg-E)  and  during  the  5-min  partition  test. 


Phenotype 

Description 

During  direct  exposure 

jumping  to  escape  from  Agg 

The  sudden  and  vigorous  leap  performed  by  the  Agg-E  mouse  as  an  apparent  attempt  to  evade  the 
Agg-initiated  attack. 

Fightback  to  Agg 

A  defensive  aggression  trait  displayed  by  the  Agg-E  mouse.  Triggered  by  the  Agg-initiated  attack,  the 
Agg-E  mouse  first  rose  up  on  its  hind  limbs  in  a  “boxing”  posture  and  swiftly  struck  out  using  a 
forepaw  toward  the  approaching  Agg.  This  appears  to  be  a  defensive-aggression  mechanism  closely 
resembling  the  fightback  phenotype  demonstrated  by  the  male  rats  stressed  by  shock  [75]. 

Freezing 

The  period  of  immobility  shown  by  Agg-E  mouse  as  the  absence  of  all  visible  movement  of  the  body 
except  the  movement  for  respiration 

During  partition  test 

Partifion  vigilance 

The  total  time  spent  per  visit  to  the  partition  zone.  The  six  inch  wide  area  in  front  of  the  subject 
animal’s  side  of  partition  is  marked  as  the  “partition  zone”.  This  parameter  may  also  imply  the  inverse 
of  the  partition  avoidance.  The  earful  mouse  is  expected  to  show  suppressed  vigilance. 

Withdrawal  to  peripheral  region 

The  time  spent  at  the  cage  peripheral  zone  in  comparison  to  the  time  spent  at  the  partition  zone 
(duration  of  time  spent  in  peripheral  zone/partition  zone).  The  peripheral  region  is  marked  by  the  six 
inch  wide  area  across  the  three  sidewalls  of  the  cage  used  for  residing  the  subject  mouse.  The  fearful 
mouse  is  expected  to  spend  more  time  in  periphery. 

Reduced  locomotion 

The  ratio  of  the  distance  traveled  inside  the  partition  zone  vs  distance  traveled  the  total  cage  portion 
available  to  the  subject  mouse  during  partition  test.  Reduced  locomotion  is  expected  to  be 
pronounced  in  stressed  animals. 

Freezing 

The  period  of  immobility  defined  by  <10%  shift  of  body  position  between  sequential  frames  of  the 
fifteen  frames  sampled  per  second.  The  ratio  of  freezing  duration  and  total  time  in  partition  zone  is 
presented.  Freezing  is  believed  to  be  positively  correlated  with  fear  response. 

Grooming 

A  behavioral  endpoint  consists  of  licking  of  paws  and  followed  by  stroking  of  nose,  head  and  body  that 
may  occur  in  any  order.  Total  grooming  duration  presented  herein  is  measured  by  visual  inspection  of 
the  videos.  It  is  believed  to  be  a  common  response  to  fear. 

Overall  activity 

The  average  velocity  of  the  mice  during  the  partition  test. 

Table  2 

The  fraction  (in  percentage)  of  the  total  Agg-E  population  showing  myocardial  degeneration  and  lymphohistiocytic  myocarditis. 


Strain 

5-day  Agg-E  stress  (%) 

10-day  Agg-E  stress  (%) 

1  d  delay 

1.5  w  delay 

1  d  delay 

4  w  delay 

C57BL/6j 

Myocardial  Degeneration 

12.5 

60.0 

36.4 

20.0 

Lymphohistiocytic  Myocarditis 

75.0 

None 

45.5 

None 

DBA/2j 

Myocardial  Degeneration 

None 

None 

None 

None 

Lymphohistiocytic  Myocarditis 

None 

None 

None 

None 

BALB/cj 

Myocardial  Degeneration 

36.4 

None 

None 

None 

Lymphohistiocytic  Myocarditis 

None 

None 

None 

None 
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Fig.  2.  Ethogram  of  Agg-E  mice  during  the  partition  test.  The  endpoints  of  Agg-E  mice  recorded  during  the  partition  test  are  normalized  against  the  age-  and  strain-matched 
controls  (Agg-E  mice/control)  and  the  mean  ±  SEM  of  Agg-E  mice  are  plotted.  The  horizontal  dotted  line  passing  through  1 .0  of  y-axis  marks  the  baseline  defined  by  the 
normalized  phenotypes  exhibited  by  the  control  mice.  Welch’s  t-test  compared  Agg-E  mice  vs.  Control:  <  0.001 ,  **p  <  0.01 ,  *p  <  0.05,  *p  <  0.1 ;  and  Welch’s  f-test  compared 

two  delay  periods  (1  d  vs.  1.5  w  for  5-day  Agg-E  stress  and  1  d  vs.  4  w  10-day  Agg-E  stress):  ***p  <  0. 001,  **p<  0.01,  *p<  0.05,  ^p<  0.1  (iV= 6-1 7  mice/group).  (A)  Partition  vigilance, 
(B)  Withdrawal  time  in  periphei'y,  (C)  Freezing,  (D)  Grooming. 
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Real-time  PCR  was  carried  out  using  the  RT^  Profiler  PCR  Array 
System  ( SA  Biosciences,  Frederick,  MD)  on  the  ABl  7900HT  platform 
(Life  Technologies)  following  the  protocol  described  earlier  [32], 

2.8.  Statistical  analysis 

Statistical  analyses  were  performed  using  GraphPad®  version 
5.0  software  (GraphPad  Software,  Inc.,  La  Jolla,  CA).  Experimental 
results  were  represented  as  mean  ±  SEM.  As  per  the  earlier  descrip¬ 
tion  [20,33],  fitting  models  were  used  to  evaluate  the  longitudinal 
trajectories  of  the  phenotypes  (such  as  body  weight,  temperature, 
and  urine  markings)  observed  over  the  5-/10-day  Agg-E  sessions. 
The  robustness  of  fit  of  the  curve  was  measured  by  and  the  p  value 
estimated  the  degree  of  deviation  of  the  fitting  curve’s  slope.  We 
limited  the  reporting  to  the  linear  regression  model  that  essentially 
explained  most  of  the  longitudinal  trajectories. 

Unpaired  t-tests  with  Welch’s  correction  were  used  to  evaluate 
the  shift  of  the  physiological  and  phenotypic  endpoints.  Pair-wise 
associations  between  phenotypic  parameters  were  evaluated  using 
the  Spearman  correlation.  Two-way  and  3-way  analysis  of  vari¬ 
ance  (ANOVA)  and  multivariate  analysis  of  variance  (MANOVA) 
were  performed.  The  principal  component  analysis  (PGA)  [34-39] 
is  reported  in  the  supplementary  data. 

We  used  R  package  version  3.0.0  (http://cran.r-project.org) 
to  compute  PCA  and  other  statistical  analysis,  and  plotted  the 
variables  using  Biplot.  The  loading  factors  contributing  most  sig¬ 
nificantly  to  the  top  three  principal  components  were  computed 
by  ANOVA  and  MANOVA.  Similar  schemes  were  used  in  the  past  to 
explain  the  ethogram  [34-39]. 

We  used  p  <  0.05  as  the  significance  threshold  for  all  the  calcu¬ 
lations. 

3.  Results 

3.1.  Strain-to-stmin  variances  in  the  longitudinal  shift  of  the 
bodyweight,  temperature  and  urine  markings  during  the  Agg-E 
stress  session 

We  measured  the  bodyweights  (Fig.  1 A  and  Fig.  S2A)  and  temper¬ 
atures  of  the  Agg-E  mice  (Fig.  S2B-C)  before  (pre-stress)  and  after 
(post-stress)  the  6  h  ‘cage-within-cage’  confinement,  and  the  urine 
markings  (Fig.  1 B)  were  measured  during  the  6  h  ‘cage-within-cage’ 
confinement.  The  pre-stress  bodyweight  and  temperature  defined 
the  baselines  of  the  respective  studies.  The  deviations  of  bodyweight 
and  temperature  from  their  respective  baselines  were  normalized 
by  the  controls  (Agg-E/control  mice)  and  the  daily  averages  of  the 
normalized  data  (±SEM)  were  plotted  (N  =  5-21  mice  per  group). 
Pair-wise  Welch’s  t-test  computed  the  differences  between  control 
and  Agg-E  mice  at  individual  time  points. 

The  linear  regression  model  was  used  to  interpret  the  longi¬ 
tudinal  trend  of  changing  the  phenotypes  as  per  earlier  reports 
[20,33,40].  C57BL/6j  showed  a  gradual  and  significant  weight  gain 
in  the  course  of  a  10-day  Agg-E  session;  in  contrast,  BALB/cj 
displayed  a  relatively  descending  trend  in  weight  significantly  dif¬ 
ferent  from  C57BL/6J  (Fig.  lA). 

A  strain-to-strain  difference  was  also  evident  in  the  longitudi¬ 
nal  shift  of  the  post-stress  bodyweights  of  Agg-E  mice  (Eig.  S2A). 
C57BL/6j  Agg-E  mice  showed  a  sharp  weight  gain  on  day  5.  Showing 
an  opposite  trend,  the  bodyweights  of  BALB/cj  Agg-E  mice  displayed 
an  early  loss  on  day  4  followed  by  a  regain  up  to  the  baseline. 
DBA/2j  mice  showed  a  significant  linear  inclination  of  post-stress 
bodyweight  over  the  10  days. 

Unlike  the  bodyweight  profile,  the  body  temperature  was  found 
to  be  unrelated  to  the  Agg-E  stress  (Fig.  S2B-C)  although  other 
studies  have  reported  the  alteration  of  core  body  temperature  as 


a  consequence  of  emotional  disturbances  [41,42].  Nevertheless,  a 
significant  strain-to-strain  difference  was  registered  in  the  longi¬ 
tudinal  alteration  of  pre-stress  temperature  of  Agg-E  mice. 

Examining  the  urine  markings  on  the  blotter  papers,  we  observed 
strain-to-strain  differences  (Eig.  1 B)  in  territorial  behavior  in  course 
of  the  1 0-day  Agg-E  sessions.  We  normalized  the  daily  urine  counts 
of  the  Agg-E  mouse  by  its  strain-specific  controls  (Agg-E/Control). 
The  numbers  of  urine  markings  of  the  Agg-E  mice  of  C57BL/6j  and 
DBA/2j  increased  with  time.  BALB/cj  Agg-E  mice,  on  the  other  hand, 
showed  a  steady  decline  of  the  territorial  behavior  with  a  terminal 
plateau.  At  day  10,  the  territorial  behavior  of  C57BL/6j  and  BALB/cj 
were  significantly  different  from  their  respective  controls  as  com¬ 
puted  by  Welch’s  t-test. 

3.2.  Strain-to-strain  variability  in  adopting  strategies  to  evade 
Agg  mice 

Eig.  IC-D  and  Fig.  S2D  depict  three  ethograms  manually 
recorded  during  the  daily  bouts  of  direct  exposures  in  the  course 
of  10-day  Agg-E  sessions.  An  individual  Agg-E  mouse  displaying 
a  particular  trait  in  at  least  two  of  three  direct  exposures  per  day 
was  classified  positive  for  that  day.  The  longitudinal  trend  of  the 
behavioral  shift  was  robustly  explained  by  the  linear  regression 
model. 

3.2.1.  Jumping  (Fig.  1C) 

A  significant  strain-to-strain  variability  was  evident  in  adopting 
jumping  as  a  mode  to  escape  from  the  Agg  mice.  In  comparison  to 
C57BL/6j  and  BALB/cj,  30%-60%  more  DBA/2j  Agg-E  mice  displayed 
jumping,  respectively.  During  the  final  two  days  of  the  10-day  Agg- 
E  session,  40%  of  DBA/2j  Agg-E  mice  displayed  a  rate  of  jumping, 
which  was  twice  that  of  the  other  two  strains. 

In  fact,  DBA/2j  Agg-E  mice  were  possibly  able  to  avoid  the  attacks 
from  the  Agg  mice  due  to  their  ability  to  jump  so  frequently.  On 
average,  Agg  mice  were  able  to  inflict  only  6  strikes/min  on  DBA/2j 
Agg-E  mice;  whereas  they  inflicted  10  strikes  on  C57BL/6J  and 
BALB/cj  Agg-E  mice  within  48  s  during  the  direct  exposures.  The 
gradual  elevation  of  the  post-stress  body  weights  of  DBA/2j  Agg-E 
mice  (Fig.  S2A)  was  oppositely  and  significantly  correlated  with  the 
diminishing  tendencies  of  jumping  (r=  -0.65). 

3.2.2.  Fightback(Fig.  ID) 

Fightback  to  the  Agg  mouse  was  another  phenotypic  endpoint 
commonly  displayed  by  all  three  strains,  though  a  significant  vari¬ 
ability  was  noted  among  the  three  strains.  This  trait  was  gradually 
and  significantly  diminished  in  C57BL/6j  and  DBA/2j ;  while  the  size 
of  BALB/cj  subpopulation  displaying  jightbacfc  remained  same.  The 
longitudinal  elevation  of  the  pre-stress  bodyweights  of  C57BL/6j 
Agg-E  mice  was  negatively  correlated  with  its  declining  fightback 
tendency  (r=-0.65).  The  fightback  tendencies’  contrasting  corre¬ 
lations  with  the  territorial  behaviors  of  DBA/2j  and  BALB/cj  were 
noteworthy.  In  DBA/2j,  the  declining  jightback  tendency  was  nega¬ 
tively  correlated  (r= -0.94)  with  its  increasing  territorial  behavior, 
whereas  in  BALB/cj  fightback  tendency  and  territorial  behavior 
showed  positive  correlation  as  both  decreased  over  time  (r=0.86). 

3.2.3.  Freezing  (Fig.  S2D) 

We  observed  significant  strain-to-strain  differences  in  freezing 
behavior.  In  fact,  the  size  of  population  displaying  freezing  signif¬ 
icantly  decreased  with  time  across  all  three  strains  of  mice.  The 
declining  tendency  of  freezing  was  positively  correlated  with  the 
decreasing  post-stress  bodyweights  of  BALB/cj  Agg-E  mice  (r= 0.72), 
but  negatively  correlated  with  the  increasing  post-stress  body- 
weights  of  DBA/2j  Agg-E  mice  (r=-0.68). 

In  the  case  of  DBA/2j  Agg-E  mice,  their  jumping,  fight¬ 
back  and  freezing  behaviors  decreased  over  time  in  significantly 
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correlated  fashions:  r  {jumping  vs.  fightback)  =  0.74-,  r  (jumping  vs. 
freezing)  =  0.90\  r  (fightback  vs.  freezing)  =  0.58.  The  jumping  and 
fightback  tendencies  of  C57BL/6J  Agg-E  mice  diminished  in  a  corre¬ 
lated  fashion  (r=0.60).  Interestingly,  no  behavioral  traits  showed 
significant  correlation  in  BALB/cj  Agg-E  mice. 

3.3.  Assessment  of  the  contributing  factors  defining  the 
behavioral  plasticity  during  the  partition  test 

We  hypothesized  that  the  post-stress  behavioral  plasticity 
was  primarily  modulated  by  the  three  following  factors:  (i)  the 
strain  differences  or  the  genotypic  heterogeneity  defined  by  Strain 
(C57BL/6j,  DBA/2j  and  BALB/cj),  (ii)  the  magnitude  of  Agg-E  stress 
defined  by  Stress  (controls  vs.  Agg-E  mice),  and  (iii)  the  duration 
of  the  delay  period  defined  by  Delay  (Id  and  1.5 w  post  5-day 
Agg-E  session,  and  1  d  and  4w  post  10-day  Agg-E  session).  The¬ 
oretically,  the  potential  roles  of  undetermined  additional  factors  in 
mediating  the  behavioral  plasticity  cannot  be  dismissed;  however 
their  impacts  should  be  Insignificant  due  to  the  carefully  controlled 
design  of  this  study.  Together,  the  Stress  and  Delay  were  defined  as 
the  Stress  history. 

Two-way  ANOVA  (Table  3  and  Table  SI  A)  captured  the  indepen¬ 
dent  and  interactive  impacts  of  strain  differences  and  the  Stress 
history.  Deconstructing  the  Stress  history  into  two  independent 
univariates,  three-way  ANOVA  (Table  3  and  Table  SIB)  captured 
the  Independent  and  interactive  contributions  of  all  three  variables 
of  Interest.  MANOVA  explained  their  interactive  contribution  by 
manipulating  the  three  as  dependent  variables  (Table  3  and  Table 
SIC). 

We  also  computed  the  impact  of  Stress  and  Delay  on  all  three 
individual  mouse  strains,  since  the  interplay  of  these  two  vari¬ 
ables  defines  the  onset  of  PTSD-like  features  [20,43],  We  previously 
reported  the  corresponding  result  manifested  only  by  C57BL/6j 
male  mice  [20]. 

Pair-wise  Welch’s  t-test  computed  the  differences  between 
control  and  Agg-E  mice  measured  at  different  time  points.  Corre¬ 
sponding  PCA  analysis  is  reported  in  the  supplementary  data. 

3.3. 1 .  Partition  vigilance  (Fig.  2A,  Table  3 ) 

5-day  Agg-E  session  vigilance:  One  day  after  the  5-day  Agg-E 
session,  Agg-E  mice  of  C57BL/6j  showed  significantly  reduced  par¬ 
tition  vigilance,  but  1.5  w  delay  significantly  increased  the  duration 
per  visit  to  the  partition.  Agg-E  mice  of  DBA/2j  and  BALB/cj  strains 
showed  control-like  vigilance  1  d  after  5-day  Agg-E  sessions.  Subse¬ 
quent  1.5w  delay  did  not  alter  the  vigilance  tendencies  of  DBA/2j 
Agg-E  mice,  but  significantly  declined  among  the  BALB/cj  Agg-E 
mice. 

The  variability  in  partition  vigilance  was  explained  by  the  inter¬ 
play  of  mouse  Strain  and  the  Stress  history.  The  exclusive  impact  of 
Stress,  the  interactive  impact  of  Strain  and  Delay,  and  the  com¬ 
prehensive  interaction  among  all  three  factors  also  emerged  as 
significant. 

MANOVA  found  no  significant  interactive  impact  of  the  three 
variables  acting  dependently. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  nearly  significant  in  C57BL/6j 
[20]  and  BALB/cj  Agg-E  mice,  but  not  in  DBA/2j  Agg-E  mice. 

10-day  Agg-E  session  vigilance:  1  d  after  the  10-day  Agg-E  session, 
the  partition  vigilance  was  significantly  suppressed  among  the 
Agg-E  mice  of  C57BL/6j  as  compared  to  controls.  The  4  w  delay  sig¬ 
nificantly  increased  the  vigilance,  returning  it  to  the  control  level. 
The  Agg-E  mice  of  DBA/2j  and  BALB/cj  displayed  control-like  par¬ 
tition  vigilance  1  d  after  the  10-day  Agg-E  sessions,  which  declined 
significantly  with  Delay. 

The  impacts  of  Strain,  Stress  history,  and  their  interplay 
were  significant.  Three-way  ANOVA  further  explained  significant 


contributions  from  Stress,  Strain  x  Stress  and  Strain  x  Delay,  and 
moderate  contributions  from  Delay  and  Its  interplay  with  Stress 
(Stress  X  Delay). 

MANOVA  analysis  suggested  significant  interactive  influence  of 
all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  Interplay 
between  Stress  and  Delay  emerged  as  highly  significant  in  the  Agg- 
E  mice  of  C57BL/6j  [20]  and  moderately  significant  in  BALB/cj  and 
DBA/2j. 

3.3.2.  Withdrawal  time  in  periphery  (Fig.  2B,  Table  3) 

5-day  Agg-E  session  withdrawal  time:  C57BL/6j  Agg-E  mice  spent 
significantly  longer  withdrawal  time  in  periphery  1  d  and  1.5  w  after 
the  5-day  Agg-E  session. 

The  withdrawal  tendency  to  the  peripheral  region  of  DBA/2j 
Agg-E  mice  was  twice  that  of  the  controls  1  d  after  the  5-day  Agg-E 
session  and  was  significantly  reduced  with  the  subsequent  1.5  w 
delay.  In  sharp  contrast,  BALB/cj  Agg-E  mice  showed  only  mini¬ 
mal  withdrawal  traits  Id  after  the  5-day  Agg-E  session,  but  the 
withdrawal  time  was  increased  significantly  after  1.5  w  delay. 

The  variance  of  this  behavioral  parameter  was  explained  by  the 
exclusive  and  Interactive  effects  of  Strain  and  Stress  history.  Three- 
way  ANOVA  found  significant  impact  from  the  combined  effect  of 
three  interactive  factors  along  with  the  exclusive  impact  of  Stress, 
and  its  interplay  with  Delay. 

MANOVA  analysis  suggested  significant  interactive  influence  of 
all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  Interplay 
between  Stress  and  Delay  emerged  as  highly  significant  in  Agg-E 
mice  of  C57BL/6j  [20],  moderately  significant  in  BALB/cj,  but  non¬ 
significant  in  DBA/2j. 

1 0-day  Agg-E  session  withdrawal  time:  One  day  after  the  1 0-day 
Agg-E  session,  C57BL/6j  Agg-E  mice  spent  significantly  longer  with¬ 
drawal  time  in  periphery;  the  trait,  however,  was  reduced  by  50% 
after  the  4w  delay.  Agg-E  mice  of  the  DBA/2j  strain  spent  consis¬ 
tently  more  time  in  the  periphery  than  the  controls  1  d  and  4  w  after 
the  1 0-day  Agg-E  session.  A  late  marginally  significant  increment  of 
withdrawal  time  to  the  periphery  was  recorded  for  BALB/cj  Agg-E 
mice. 

The  exclusive  impact  of  the  Stress  history  and  its  interplay  with 
Strain  emerged  significant.  The  contribution  from  Stress  history 
was  primarily  explained  by  Stress.  The  interplay  between  Strain 
and  Delay  and  the  combined  interplay  among  all  three  variables 
was  significant. 

MANOVA  analysis  suggested  that  there  was  significant  interac¬ 
tive  influence  of  all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  Interplay 
between  Stress  and  Delay  emerged  as  highly  significant  in  Agg-E 
mice  of  C57BL/6j  [20],  but  non-significant  in  BALB/cj  and  DBA/2j. 

3.3.3.  Freezing  (Fig.  2C,  Table  3 ) 

5-day  Agg-E  session  freezing:  The  Agg-E  mice  of  C57BL/6j  froze 
more  frequently  than  the  control  1  d  after  the  5-day  Agg-E  session, 
but  displayed  only  control-llke /reezing  1.5  w  post  5-day  Agg-E 
session.  The  Agg-E  mice  of  DBA/2j  and  BALB/cj  displayed  control¬ 
like  freezing  1  d  after  5-day  Agg-E  session.  Delay,  however,  caused 
a  significant  reduction  of  freezing. 

Only  the  exclusive  contribution  of  Strain  emerged  as  significant. 
MANOVA  found  non-significant  contribution. 

Considering  the  three  strains  individually,  the  Interplay 
between  Stress  and  Delay  emerged  as  non-slgnlficant  across  the 
strains. 

10-day  Agg-E  session  freezing:  The  Agg-E  mice  of  C57BL/6j 
showed  two-fold  higher  freezing  than  the  controls  1  d  after  the 
1 0-day  Agg-E  session;  and  4  w  delay  diminished  their  freezing  dura¬ 
tions  to  control  levels.  The  freezing  durations  of  Agg-E  mice  of 


Table  3 

ANOVA  analyses  of  the  behavioral  parameters  investigated  during  the  partition  test. 


Factor 

Partition  vigilance 

Withdrawal  time  in  periphery 

Freezing 

Grooming 

Reduced  locomotion 

Overall  activity 

p  value 

F  value 

p  value 

F  value 

p  value 

F  value 

p  value 

F  value 

p  value 

F  value 

p  value  F  value 

2-way  ANOVA 

5-day  Agg-E  stress  followed  by  1  d  or  1 .5  w  delay 

Strain 

- 

- 

<0.0001 

25.03 

0.04 

3.23 

- 

- 

- 

- 

- 

Stress  history 

- 

- 

<0.0001 

16.91 

- 

- 

- 

- 

0.001 

3.93 

- 

Strain  x  Stress  history 

0.02 

2.57 

<0.0001 

7.09 

- 

- 

- 

- 

- 

- 

0.07  1.96 

10-day  Agg-E  stress  followed  by  1  d  or  4  w  delay 

Strain 

0.02 

3.88 

- 

- 

<0.0001 

11.70 

- 

- 

- 

- 

- 

Stress  histoi'y 

0.002 

5.07 

0.0005 

6.37 

0.0001 

7.17 

0.0008 

5.86 

0.0001 

7.27 

- 

Strain  x  Stress  history 

0.0003 

4.53 

0.03 

2.42 

<0.0001 

7.39 

0.01 

2.88 

- 

- 

- 

3-way  ANOVA 

5-day  Agg-E  stress  followed  by  1  d  or  1 .5  w  delay 


Strain 

- 

- 

<0.0001 

25.03 

Stress 

0.04 

4.28 

<0.0001 

22.16 

Delay 

- 

- 

- 

- 

Strain  x  Stress 

- 

- 

- 

- 

Strain  x  Delay 

0.05 

2.99 

0.04 

3.42 

Stress  X  Delay 

- 

- 

- 

- 

Strain  x  Stress  x  Delay 

0.03 

3.48 

0.03 

3.46 

10-day  Agg-E  stress  followed  by  1  d  or  4  w  delay 
Strain 

0.01 

11.71 

Stress 

0.004 

8.53 

<0.0001 

20.17 

Delay 

0.08 

3.17 

- 

- 

Strain  x  Delay 

0.01 

4.52 

- 

- 

Strain  x  Delay 

0.005 

5.33 

0.01 

4.38 

Stress  X  Delay 

0.09 

2.90 

- 

- 

Strain  x  Stress  x  Delay 

0.006 

5.18 

0.02 

4.24 

0.04 

3.23 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.001 

11.27 

0.03 

5.04 

0.05 

2.97 

- 

- 

- 

- 

0.09 

2.47 

<0.0001 

11.70 

<0.0001 

21.47 

0.0007 

11.86 

<0.0001 

21.79 

- 

- 

<0.0001 

10.29 

_ 

_ 

_ 

_ 

_ 

_ 

0.003 

6.22 

0.0004 

8.34 

- 

- 

- 

- 

- 

- 

0.02 

5.77 

- 

- 

- 

- 

0.004 

5.74 

0.07 

2.65 

- 

- 

- 

- 
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DBA/2j  did  not  change  from  controls  1  d  and  4  w  after  10-day  Agg- 
E  session.  BALB/cj  Agg-E  mice  displayed  a  modestly  long  freezing 
duration  1  d  after  10-day  Agg-E  session,  which  was  increased  sig¬ 
nificantly  after  the  4  w  delay. 

The  exclusive  and  interactive  impacts  of  Strain  and  Stress  history 
emerged  as  significant.  The  interplay  of  Strain  with  Stress  and  Delay 
respectively  emerged  as  significant:  in  addition,  the  combinatory 
influence  of  the  three  variables  was  significant. 

MANOVA  analysis  suggested  significant  interactive  influence  of 
all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  highly  significant  in  Agg-E 
mice  of  C57BL/6j  [20]  and  BALB/cj;  but  non-significant  in  DBA/2J. 

3.3.4.  Grooming  (Fig.  JD,  Table  3) 

5-day  Agg-E  session  grooming:  Among  three  mouse  strains,  only 
the  Agg-E  mice  of  C57BL/6J  showed  a  delay-induced  significant  shift 
in  grooming  duration. 

None  of  variables  was  significant  in  interpreting  grooming 
as  inferred  from  2-/3-way  ANOVA.  MANOVA  analysis,  however, 
suggested  modest  interactive  influence  of  all  three  dependent  vari¬ 
ables. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  non-significant  in  modulat¬ 
ing  grooming  across  the  strains. 

1 0-day  Agg-E  session  grooming:  Grooming  times  of  Agg-E  mice  of 
both  C57BL/6j  and  DBA/2j  strains  were  nearly  significantly  higher 
than  the  controls  1  d  and  4w  post  10-day  Agg-E  session.  On  the 
other  hand,  the  Agg-E  mice  of  BALB/cj  strain  spend  significantly 
less  time  grooming  than  the  controls  1  d  after  the  stress.  After 
4w  delay,  however,  the  grooming  time  significantly  increased  in 
BALB/cj. 

The  exclusive  impact  of  Stress  history  and  its  interplay  with 
Strain  emerged  as  significant.  The  contribution  from  Stress  bistory 
was  primarily  explained  by  Stress.  The  interplay  between  Strain 
and  Delay  and  the  combinatory  influence  of  all  three  factors  were 
significant. 

MANOVA  analysis  suggested  significant  interactive  influence  of 
all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  highly  significant  in  BALB/cj 
and  non-significant  in  C57BL/6j  and  DBA/2j. 

3.3.5.  Reduced  locomotion  (Eig.  S3A,  Table  3) 

5-day  Agg-E  session  locomotion:  The  Agg-E  of  C57BL/6j  showed 
significantly  decreased  locomotion  1  d  after  5-day  Agg-E  session. 
The  Agg-E  mice  of  DBA/2j  displayed  control-like  locomotion  both 
1  d  and  1.5  w  after  the  5-day  Agg-E  session.  The  Agg-E  mice  of  the 
BALB/cj  strain,  on  the  other  hand,  while  normal  at  1  d,  displayed 
significantly  less  locomotion  after  1.5  w  delay.  Thus,  the  stress- 
induced  inhibition  of  locomotion  was  delayed  in  BALB/cj,  compared 
to  the  other  two  strains. 

MANOVA  analysis  found  significant  interactive  contribution 
among  all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  highly  significant  in  BALB/cj 
Agg-E  mice,  but  non-significant  in  C57BL/6j  and  DBA/2j. 

1 0-day  Agg-E  session  locomotion:  The  locomotion  of  Agg-E  mice 
of  C57BL/6j  strain  showed  significant  diminution  1  d  after  the  10- 
day  Agg-E  session;  the  4w  delay  caused  a  modest  recovery.  The 
Agg-E  mice  of  the  DBA/2j  strain  displayed  decreased  locomotion 
both  1  d  and  4w  after  the  10-day  Agg-E  session,  compared  to  the 
controls.  Moreover,  the  Agg-E  mice  of  the  DBA/2j  strain  displayed 
further  diminution  of  locomotion  after  the  4w  delay.  The  Agg-E 
mice  of  the  BALB/cj  strain  consistently  showed  more  decreased 


locomotion  compared  to  controls  both  1  d  and  4  w  after  the  10-day 
Agg-E  session. 

Delay  post  stress  was  the  significant  factor  for  interpreting  the 
variability  in  decreased  locomotion. 

MANOVA  analysis  suggested  significant  interactive  influence  of 
all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  non-significant  across  tbe 
strains. 

3.3.6.  Overall  activity  (Fig.  S2B,  Table  3) 

5-day  Agg-E  session  Overall  activity:  The  overall  activities  of 
C57Bl/6j  Agg-E  mice  were  similar  to  the  controls  1  d  after  the  5-day 
Agg-E  session,  and  decreased  marginally  after  1 0  d  delay.  In  DBA/2j, 
the  activity  level  was  significantly  lower  than  the  controls  1  d  after 
the  5-day  Agg-E  session,  and  subsequent  1 .5  w  delay  resulted  in  tbe 
activity  level  returned  to  tbe  baseline.  As  in  the  case  of  C57BL/6j, 
Agg-E  mice  of  BALB/cj  displayed  control-like  activities  1  d  after  the 
5-day  Agg-E  session,  but  BALB/cj  showed  a  significant  decrease  in 
activity  than  controls  after  1.5  w  delay. 

Stress  history  particularly  defined  by  Stress  and  its  interplay 
with  Strain  displayed  a  modest  effect  on  overall  activity. 

MANOVA  analysis  suggested  moderate  interactive  influence  of 
all  three  dependent  variables. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  modestly  significant  in 
C57BL/6j  [20]  and  non-significant  in  the  other  two  strains. 

10-day  Agg-E  session  Overall  activity:  The  Agg-E  mice  of  both 
C57BL/6j  and  DBA/2j  displayed  control-like  overall  activities  1  d 
as  well  as  4w  after  the  10-day  Agg-E  sessions.  The  Agg-E  mice 
of  BALB/cj  strain  displayed  variable  but  not  significantly  elevated 
activities  compared  to  controls  1  d  after  the  10-day  Agg-E  sched¬ 
ule  but  the  increased  activity  of  BALB/cj  mice  compared  to  controls 
reached  statistical  significance  after  4  w  delay. 

None  of  the  factors  and  their  interplays  was  significant  in 
interpreting  overall  activity  as  inferred  by  2-/3-way  ANOVA  and 
MANOVA. 

Considering  the  three  strains  individually,  the  interplay 
between  Stress  and  Delay  emerged  as  non-significant  across  tbe 
strains. 

3.4.  Cardio-histopathological  analysis  of  the  Agg-E  mice 

We  examined  the  heart  tissue  for  two  potential  indicators 
of  cardiac  dysfunctions  typically  induced  by  stress,  myocar¬ 
dial  degeneration  and  lymphohistiocytic  myocarditis  (N  =  5  per 
group)  [21  [.  More  than  12%  of  the  C57BL/6j  Agg-E  mice  dis¬ 
played  myocardial  degeneration  1  d  after  the  5-day  Agg-E  session. 
This  myopathy  became  more  pronounced  after  1.5  w  delay  as 
60%  of  the  Agg-E  mice  sustained  myocardial  degeneration.  10- 
day  Agg-E  session  caused  myocardial  degeneration  among  36.4% 
of  C57BL/6j  Agg-E  mice.  Four  week  delay  may  have  permit¬ 
ted  some  reversal  of  the  atrophy,  as  only  20%  of  the  Agg-E 
mice  displayed  myocardial  degeneration.  We  diagnosed  lym¬ 
phohistiocytic  myocarditis  in  75%  and  45.5%  of  C57BL/6j  Agg-E 
mice  Id  after  the  5-day  and  10-day  Agg-E  sessions,  respec¬ 
tively.  The  subsequent  delays  might  have  positive  impact  on  this 
inflammation. 

The  cardiac  health  of  DBA/2j  mice  was  found  least  susceptible 
to  this  type  of  stress  as  none  of  the  DBA/2j  Agg-E  mice  displayed 
any  cardiac  dysfunction  after  the  Agg-E  schedule. 

One  day  after  the  5-day  Agg-E  session,  36.4%  of  BALB/cj  Agg-E 
mice  displayed  myocardial  degeneration,  but  none  of  them  showed 
lymphohistiocytic  myocarditis.  No  cardiac  myopathy  was  diag¬ 
nosed  after  the  delay  periods. 
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Fig.  3.  qPCR  result  of  the  selected  panel  of  genes.  The  fold  changes  (Agg-E 
mice/Control)  of  six  genes  are  reported  in  log2  scale.  The  horizontal  dotted  line 
passing  through  1.0  of  Y-axis  marks  the  cut-off  at  fold  change  ±1.5.  Welch’s  correc¬ 
tion  compared  Agg-E  mice  vs.  control  mice:  *p<0.05;  1  d  vs.  4w:  @  0.05<p<0.1, 
•p  <  0.05,  "p  <  0.01 ,  -"p  <  0.0001 . 

3.5.  Real-time  PCR  results  to  characterize  the  selective  genetic 
responses  to  repeated  Agg-E  stress 

We  measured  the  regulations  of  six  genes  using  qPCR,  namely 
Bdnf,  Ngf, Zwint,  Cckbr,  Slc6a4  and  FkbpS  (Fig.  3 ). Table  4  documents  a 
brief  literature  survey  linking  these  molecules  to  PTSD  and  synaptic 
plasticity.  For  this  study,  we  used  hemibrains  collected  from  Agg-E 
(N  =  4-8)  and  control  mice  (N  =  4-10)  of  all  three  strains  subjected  to 
1 0-day  Agg-E  sessions  and  euthanized  either  1  d  or  4  w  post-stress. 

Bdnf:  One  day  after  a  1 0-day  Agg-E  session,  Bdnf  was  about  four¬ 
fold  suppressed  in  C57BL/6J  Agg-E  mice.  Eour  weeks  delay  shifted 
this  to  nearly  two-fold  down-regulated  from  the  control  level.  In 


DBA/2J,  Bdnf  was  elevated  by  nearly  two-fold  among  Agg-E  mice 
1  d  after  the  10-day  Agg-E  session  and  remained  at  this  level  4w 
post-stress.  We  observed  a  2.5-fold  elevation  of  Bdnf  in  BALB/cJ 
Agg-E  mice  1  d  after  the  10-day  stress  that  reduced  significantly 
4w  post-stress. 

Ngf:  In  C57BL/6J  Agg-E  mice,  the  down-regulated  expression  of 
Ngf  1  d  after  the  10-day  Agg-E  session  became  up  regulated  after 
4w  delay.  Two-fold  elevated  expressions  were  observed  in  Agg-E 
mice  of  the  other  two  strains  1  d  post-Agg-E  session  that  declined 
to  the  control  level  after  4  w  delay. 

Zwint:  In  C57BL/6j  Agg-E  mice,  Zwint  was  significantly  sup¬ 
pressed  Id  after  the  10-day  Agg-E  session  and  was  elevated 
to  control  level  4w  post-stress.  A  control-like  expression  was 
observed  in  DBA/2j  Agg-E  mice.  A  consistently  and  significantly 
elevated  expression  of  Zwint  was  observed  among  BALB/cj  Agg-E 
mice. 

Cckbr:  One  day  after  the  1 0-day  Agg-E  session,  Cckbr  was  almost 
significantly  elevated  in  C57BL/6J  Agg-E  mice,  and  after  4w  it 
remained  2.5-fold  elevated  from  the  control.  A  late  and  nearly  two¬ 
fold  elevation  was  observed  in  DBA/2J  Agg-E  4w  after  the  10-day 
Agg-E  session.  There  was  a  consistent  two-fold  elevation  in  BALB/cj 
Agg-E  mice  1  d  and  4  w  after  the  10-day  Agg-E  session. 

Slc6a4:  Consistent  suppression  of  Slc6a4  was  observed  in 
C57BL/6j  Agg-E  mice,  and  even  more  so  after  4w  delay.  In  con¬ 
trast,  a  control-like  expression  was  observed  in  DBA/2j  Agg-E  mice 
1  d  and  4  w  after  the  10-day  Agg-E  session.  Among  BALB/cj  Agg-E 
mice,  we  found  control-like  regulation  1  d  after  the  10-day  Agg-E 
session  followed  by  slight  suppression  4  w  later. 

FkbpS:  A  consistent  two-fold  suppression  of  FkbpS  was  regis¬ 
tered  Id  and  4w  post-lO-day  Agg-E  session  in  C57BL/6j  Agg-E 
mice.  In  DBA/2j  Agg-E  mice,  FkbpS  remained  two-fold  elevated 
1  d  and  4w  post-1 0-day  Agg-E  session.  A  slight  suppression  was 
observed  in  BALB/cj  Agg-E  mice  4  w  post-1 0-day  Agg-E  session. 


Table  4 

A  brief  literature  survey  linking  the  panel  of  genes  of  interest  to  PTSD  and  co¬ 
morbidities. 


Brain  derived  neurotrophic 
factor  (BDNF) 


Nerve  growth  factor  (NGF) 


ZWIO  interactor  (ZWINT) 


Cholecystokinin  B  receptor 
(CCKBR) 


Solute  carrier  family  6 
(neurotransmitter 
transporter,  serotonin), 
member  4  (SLC6A4) 

FK506  binding  protein  5 
(FKBPS) 


Human:  Recent  PTSD  patients  have  higher 
BDNF  than  long  term  PTSD  patients  likely  due 
to  the  biased  consolidation  of  fear  memory 
[55]. 

Rat:  Contextual  fear  causes  exon-specific  BDNF 
expression  [66]. 

Rat:  Sleep  deprivation,  western  diet  and 
predator  scent  stress  suppress  hippocampal 
BDNF  [62,63]. 

Mouse;  NGF  depletion  contributes  to  age 
related  cognitive  dysfunction  [67]. 

Mouse:  NGF  in  hippocampus  and  amygdala  is 
negatively  related  to  anxiety  but  positively 
related  to  the  conditioned  fear  [68]. 
Mammalian  cells:  ZWINT  is  involved  in 
pre-synaptic  events  in  vitro  [69]. 

Rat:  ZWINT  helps  in  neuroprotection  against 
sub-lethal  stress  in  vitro  [70]. 

Fluman:  CCKBR  is  associated  with  panic 
disorder  [59]. 

Mouse:  CCKBR  is  responsible  for 
anxiety-related  behavior  [76,77]. 

Human:  Suppressed  polymorphism  in  the 
serotonin  transporter  gene  is  associated  with 
PTSD  [58]. 

Mouse;  Suppressed  serotonin  transporter  gene 
elevates  social  avoidance  and  depression  [78]. 
Human:  Interaction  of  FKBPS  with  acute  child 
abuse  is  a  putative  predictor  of  adult  PTSD  [56]. 
Human:  Four  single  nucleotide  polymorphisms 
located  within  the  FKBPS  are  associated  with 
PTSD  risk  [57]. 

Mouse;  FKBP5  is  associated  with  chronic  stress 


172], 


4.  Discussion 

Previous  studies  exploited  genetic  variances  among  the  var¬ 
ious  mouse  strains  to  understand  its  contributions  in  shaping 
the  social  behaviors  [44],  accomplishing  the  behavioral  tasks  [45] 
and  responding  to  various  types  of  stress,  such  as  brief  foot- 
shock  [14,28],  predator’s  scent  [29]  and  early  life  trauma  [46]. 
The  dual  effects  of  the  genetic  heterogeneity  and  the  phenotypic 
challenges  (such  as  changing  stressor  from  neurological  injury 
to  emotional  assault,  introducing  environmental  novelty)  on  the 
rodents’  response  profile  (such  as  fear  plasticity,  neurotransmis¬ 
sion)  were  clearly  evident.  In  this  context,  the  present  objective 
is  to  characterize  the  interplay  between  the  genetic  heterogeneity 
and  environmental  variables  in  developing  PTSD-like  features  in 
murine  model. 

We  used  a  conspecific  Agg-E  stress  model,  a  modified  version 
of  the  aggressor-intruder  social  stress  model  [20]  potentially  sim¬ 
ulating  the  unique  experiences  of  the  soldiers,  who  repeatedly 
encounter  the  life  threatening  trauma  in  combat  theater.  Carefully 
designing  the  conspecific  and  repeated  interaction  with  physically 
threatening  traumas,  this  model  may  replicate  some  aspects  of 
combat  stress  potentially  more  effectively  than  the  conventional 
PTSD  models  such  as  the  brief  footshock  model  (not  a  prolonged 
and  repeated  stress  model)  and  predator  scent  model  (not  a  con¬ 
specific  model).  The  concern  [12,26]  that  repeated  social  stress 
may  elicit  habituation  rather  than  triggering  the  PTSD-like  features 
was  mitigated  by  carefully  implementing  the  ‘randomness’  in  the 
occurrence  of  the  life-threatening  events.  Thereby  the  ’ethological 
validity’  of  the  model  namely  the  unpredictable  and  uncontrollable 
nature  of  the  PTSD  simulating  trauma  was  established.  Moreover, 
the  ‘dose-response’-type  feature  was  presented  as  10-day  long 
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stress  in  comparison  to  5-day  long  stress  emerged  more  proficient 
in  eliciting  PTSD-like  features  [20-23],  In  addition,  the  ‘face  valid¬ 
ity’  of  PTSD-like  features  [12,26]  was  manifested  by  introducing 
the  stressed  mice  to  the  contextual  cues  immediately  after  with¬ 
drawing  the  trauma  and  at  a  delayed  interval,  wherein  4  weeks  is 
arguably  equivalent  to  nearly  3  years  of  human  life.  Supporting  our 
claim  we  demonstrated  many  traumatic  features  emerging  after 
a  delay  following  the  withdrawal  of  stress  [20-23],  Concluding, 
the  present  data  display  many  PTSD-simulatlng  phenotypic  and 
molecular  signatures  immediately  after  the  trauma-withdrawal; 
and  many  of  these  traits  were  sustained  through,  or  even  emerged 
after  the  interval  of  delay.  As  a  ‘construct  validity’  of  the  PTSD  rodent 
model  many  of  the  reported  putative  genomic  [21],  transcripto- 
mic  [22],  metabolomic  [23]  and  cellular  [20]  signatures  of  PTSD 
were  found  perturbed  in  the  present  Agg-E  repeated  stress  model. 
Accordingly,  this  rodent  model  meaningfully  simulates  the  combat 
related  PTSD  paradigm  and  validated  some  of  the  essential  criteria 
of  PTSD  rodent  models. 

Employing  this  model,  we  explored  genetic  heterogeneity  of 
three  in-bred  strains  of  mice  for  stress-elicited  display  of  PTSD- 
like  phenotypes.  Sub-population  analysis  based  on  the  Individual 
variability  in  responding  to  Agg-E  stress  is  however  beyond  the 
scope  of  present  study.  The  sample  size  used  herein  was  optimum 
for  the  phenotypic  assays,  however  the  size  of  the  subgroups  were 
often  below  to  draw  any  biologically  meaningful  Inferences;  there¬ 
fore  carrying  out  such  computation  was  overruled  for  the  present 
purpose. 

4.3.  DBAflj  displayed  a  different  strategy  to  cope  with  stress 

Compared  to  C57BL/6j,  DBA/2j  displayed  an  alteration  of  pheno¬ 
types,  such  as  jumping,  fightback  and  freezing  with  of  10-day  Agg-E 
stress.  Consistently  robust  jumping  helped  DBA/2J  to  evade  Agg 
attack,  which  was  accompanied  with  a  gradual  decrease  of  fight- 
back,  to  display  a  very  different  strategy  to  evade  Agg  mice. 

Bodyweight  was  a  superior  physiological  signature  of  Agg-E 
stress  compared  to  temperature.  The  pre-stress  bodyweight  is  the 
preferred  signature  probably  for  the  following  reasons.  The  plas¬ 
ticity  of  the  post-stress  bodyweight  could  be  likely  due  to  the 
immediate  effect  of  6  h  starvation  than  potential  metabolic  irregu¬ 
larities  caused  by  any  adjustment  of  the  rodent’s  dietary  practice. 
C57BL/6j  and  BALB/cj  mice  displayed  opposite  trends  in  both 
pre-and  post-stress  bodyweights,  supporting  a  past  report  [25].  Sig¬ 
nificant  strain-to-strain  variability  in  bodyweight  accounted  for  the 
impact  of  genetic  predisposition  in  controlling  metabolism  under 
stress. 

Decreasing  territorial  behavior  correlated  with  waning  fightback 
could  be  due  to  the  nature  of  BALB/cj  Agg-E  mice  conceding  to  the 
aggressor.  In  contrast,  DBA/2j’s  display  of  increasing  territoriality 
might  have  signified  its  improved  assessment  of  the  threat;  and 
simultaneous  reduction  of  fightback  suggested  that  It  potentially 
adopted  a  different  strategy  to  evade  the  Agg  mice  compared  to 
the  other  strains. 

4.2.  C57BL/6j  displayed  pronounced  cardiomyopathy 

A  previous  study  investigating  inter-strain  cardiac  health  found 
C57BL/6j  in  particular  was  an  outlier  in  demonstrating  significant 
susceptibility  to  hypoxic  challenges  [47].  The  causal  relation¬ 
ship  between  hypoxia  and  depression,  a  major  PTSD-comorbidity, 
should  be  noted  in  this  context  [48],  In  accordance  with  these  find¬ 
ings,  we  observed  profound  adverse  impacts  of  the  Agg-E  stress  on 
the  C57BL/6j  Agg-E  mice,  and  its  genomic  consequences  have  been 
investigated  by  us  in  a  previous  study,  limited  to  just  one  strain  [21  ]. 
Five-day  Agg-E  sessions  caused  a  late  augmentation  of  myocar¬ 
dial  degeneration,  while  10-day  Agg-E  sessions  caused  a  rather 


consistent  cardiomyopathy.  The  increased  bodyweight  of  C57BL/6j 
could  have  a  causal  relationship  with  its  poor  cardiac  health, 
as  previous  studies  suggested  close  associations  among  obe¬ 
sity/metabolic  disorder,  cardiac  health  and  PTSD  symptoms  [49,50], 
In  the  present  study  we  observed  75%  of  C57BL/6j  mice  devel¬ 
oped  lymphohistiocytic  myocarditis,  and  none  of  DBA/2j  did  so, 
and  while  36.4%  of  BALB/cj  showed  an  Immediate  effect  after  10  d 
Agg-E  stress. 


4.3.  Three  strains  displayed  different  patterns  of  coping  with 
Agg-E  stress  over  time 

Fixation  to  a  past  traumatic  event  and  reemergence  of  anx¬ 
iety,  depression  and  other  stress-associated  phenotypes  elicited 
by  contextual  provocation  are  some  of  the  major  traits  of  PTSD 
[12],  which  were  simulated  In  the  present  study  by  monitor¬ 
ing  six  behavioral  endpoints  (Table  1)  during  the  partition  test 
[20],  The  partition  vigilance  and  withdrawal  time  in  the  periphery 
counterbalances  the  curiosity  against  the  tendency  of  dissoci¬ 
ation  from  potential  threat.  The  variances  in  both  traits  after 
5-day  and  10-day  Agg-E  sessions  were  explained  by  the  significant 
interaction  between  tbe  strain  differences  and  the  Stress  his¬ 
tory.  A  delayed  Increase  in  partition  avoidance  observed  in  BALB/cj 
potentially  captured  tbe  time-evolved  emergence  of  the  rodent’s 
fear  response  to  the  traumatic  cue.  Since  the  delayed  emergence 
(or  re-experlenclng)  of  trauma  is  an  essential  criterion  of  PTSD 
[51,52],  the  delayed  onset  of  PTSD-like  features  are  of  particular 
interest. 

Grooming,  an  accepted  signature  of  anxiety,  depression  and 
hyperarousal  [53]  and  was  altered  by  the  Interaction  between 
Strain  and  Stress  history  after  the  10-day  Agg-E  session.  Re¬ 
experiencing  the  contextual  fear,  BALB/cj  Agg-E  mice  spent 
significantly  more  time  grooming  4w  after  the  10-day  Agg-E 
session. 

Freezing  and  reduced  locomotion  are  the  dual  acts  of  fear 
response,  possibly  caused  by  an  overcautious  approach  to  poten¬ 
tial  danger  and  the  tendency  of  conceding  to  the  traumatic  context 
[54],  The  significant  interplay  between  strain  difference  and  Stress 
history  explained  the  variance  in  freezing  after  the  10-day  Agg-E 
session.  Once  again,  BALB/cj  Agg-E  mice  spent  significantly  more 
time  freezing  4  w  after  the  10-day  Agg-E  session. 

It  is  important  to  note  that  the  overall  activities  of  Agg-E  mice 
were  almost  control-llke  except  for  BALB/cj  mice  whose  possible 
‘anxious’  behavior  manifested  by  frequent  grooming  might  have 
contributed  to  their  apparent  hyperactivity.  The  control-like  activ¬ 
ities  displayed  by  the  other  two  strains  in  general  suggested  that 
their  hypoactivities,  such  as  freezing  and  retarded  locomotion,  were 
potentially  driven  by  the  contextual  fear,  not  by  lethargy  or  habit¬ 
uation. 

Mapping  of  these  behavioral  parameters  In  PCA  (Fig.  S4  and 
Tables  S2-3)  revealed  similar  longitudinal  trajectories  of  coping 
with  the  stress  manifested  by  C57BL/6j  and  DBA/2j  post-5-day 
Agg-E  session,  whereas,  BALB/cj  displayed  a  reversed  longitudinal 
trajectory.  DBA/2j  Agg-E  mice  maintained  a  minimal  behavioral 
shift  post-1 0-day  Agg-E  session,  depicting  a  delay-independent 
response  to  Agg-E  stress. 

The  behavioral  changes  between  day  1  and  day  10  were  stri¬ 
kingly  different  between  C57BL/6j  and  BALB/cj  for  both  vigilance 
and  freezing  -  almost  mirror  images  of  one  another.  Vigilance  was 
increasing  in  C57BL/6j  but  decreasing  in  BALB/cj.  In  marked  con¬ 
trast,  freezing  was  decreasing  in  C57BL/6j  as  it  was  increasing  in 
BALB/cj.  The  reemergence  or  delayed  development  of  many  psychi¬ 
atric  symptoms  is  typical  among  PTSD  patients  [51,52],  therefore 
these  distinctive  longitudinal  patterns  of  behavioral  shifts  demon¬ 
strated  by  BALB/cj  would  be  important  to  examine  further. 
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4.4.  C57Bl/6j  and  BALB/cj  displayed  perturbed  expression  of 
selected  PTSD  signatures 

Genes  associated  with  synaptic  plasticity  {Bdnf,  Ngf  and  Zwint) 
and  genes  encoding  receptors  for  cholcystokinin  (e.g.  Cckbr),  co¬ 
chaperone  of  the  glucocorticoid  (e.g.  FkbpS)  and  transporter  of 
serotonin  (e.g.  Slc6a4)  were  evaluated.  Potential  association  of 
Bdnf  [55],  FkbpS  [56,57]  and  Slc6a4  [11,58]  with  PTSD  has  been 
established,  while  Cckbr  has  been  linked  to  anxiety  [59],  The  mea¬ 
surements  in  the  present  study  were  limited  for  being  focused  on 
the  hemibrains,  not  on  the  different  brain  regions.  Furthermore,  the 
outcome  could  be  potentially  compromised  by  the  collection  of  the 
small  midline  structures  such  as  thalamic  and  hypothalamic  nuclei. 

The  complex  relationship  of  Bdnf  with  synaptic  plasticity, 
depression  and  memory  consolidation  [60,61  ]  might  be  taken  into 
account  to  explain  its  variable  expression  across  the  strains.  We 
observed  suppressed  Bdnf  in  C57BL/6j  in  accordance  with  the  stud¬ 
ies  probing  sleep  deprivation  [62],  predator  stress  [63],  acute  social 
defeat  [64]  and  social  isolation  [65],  Exon-specific  elevation  of  Bdnf 
in  rat  hippocampus,  however,  was  reported  as  a  consequence  of 
restraint  stress  [66],  In  comparison  to  patients  with  long-standing 
PTSD,  patients  with  recently  acquired  PTSD  have  been  diagnosed 
with  bigber  serum  Bdnf  which  was  interpreted  as  a  causal  con¬ 
sequence  of  the  biased  consolidation  of  fear  memory  [55],  This 
observation  may  explain  the  elevated  Bdn/found  immediately  after 
10-day  Agg-E  stress  in  BALB/cj,  followed  by  delayed  suppression. 
Hence,  frequent  manifestation  of  delayed  arousal  of  PTSD-like  traits 
in  BALB/cj  such  as  withdrawal  to  periphery,  freezing,  and  groom¬ 
ing  could  be  attributed  to  the  biased  consolidation  of  fear  memory 
causally  related  to  the  suppression  of  Bdnf  in  hemibrain. 

Significant  changes  of  Ngf  were  limited  to  C57BL/6j,  where 
there  was  a  decrease  on  day  1  followed  by  an  increase  on  day 
10.  Initial  suppression  of  Ngf  in  C57BL/6j  might  hypothetically  be 
associated  with  cognitive  dysregulation  [67],  increased  anxiety  and 
poor  fear  conditioning  [68]  as  the  immediate  effect  of  10-day  Agg- 
E  stress.  Subsequent  elevation  of  Ngf  after  4w  delay,  associated 
with  decreased  withdrawal  and  decreased  freezing  putatively  sug¬ 
gests  the  molecule’s  association  with  reversal  of  fear  conditioning 
[68],  Zwint,  another  gene  closely  associated  with  synaptic  plas¬ 
ticity  [69,70],  was  suppressed  but  only  on  day  one  in  C57BL/6j. 
Zwint  showed  a  significant  elevation  at  both  time  points  in  BALB/cj, 
indicating  a  possible  strain-specific  and  stress-induced  prolonged 
alteration  of  synaptic  plasticity.  In  accordance  with  their  endophe- 
notypes,  DBA/2j  displayed  minimum  perturbation  of  Bdnf,  Ngf  and 
Zwint,  which  might  be  associated  with  its  apparent  resilience  in 
the  paradigm.  Onset  of  stress  typically  diminishes  the  activity  of 
the  co-chaperone  domain  F1<BP5  [71],  which  has  been  shown  in 
the  chronic  stress  model  [72]  and  presented  in  the  Agg-E  model 
assessing  C57BL/6j  and  BALB/cj. 

4.5.  Summary 

The  patterns  of  response  to  Agg-E  stress  emerged  significantly 
distinct  among  the  strains.  C57Bl/6j  showed  very  poor  cardiac 
health  and  a  very  different  phenotypic  strategy  of  coping  with  Agg- 
E  stress,  compared  to  other  strains.  Delayed  emergence  of  many 
PTSD-simulating  features  distinguished  BALB/cj  from  the  other 
two  strains,  with  exception  of  myocardial  vulnerability,  which  was 
expressed  early.  The  strain-to-strain  variability  in  characterizing 
the  PTSD-like  phenotypes  was  primarily  explained  by  the  inter¬ 
play  of  the  genetic  heterogeneity  interacting  with  stress  history. 
The  contribution  of  Stress  history  was  primarily  defined  by  the 
magnitude  of  the  stress  and  its  interplay  with  the  delay  period 
(Stress  X  Delay). 

It  is  further  interesting  to  note  that  there  is  a  close  associa¬ 
tion  of  the  regulation  of  a  selected  pool  of  genes  with  phenotypic 


plasticity.  Specifically,  the  comparatively  distant  molecular  expres¬ 
sions  seen  in  DBA/2j  were  in  contrast  with  the  other  two 
strains.  Studies  identifying  phenotype-specific  molecular  foot¬ 
prints  [66,73,74]  are  underway:  more  comprehensive  approaches, 
with  a  focus  on  investigating  individual  brain  regions  in  respect  to 
PTSD  models  are  warranted. 

Disclaimer 

The  views,  opinions,  and/or  findings  contained  in  this  report  are 
those  of  the  authors  and  should  not  be  construed  as  official  Depart¬ 
ment  of  the  Army  position,  policy,  or  decision,  unless  so  designated 
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Research  was  conducted  in  compliance  with  the  Animal  Wel¬ 
fare  Act,  and  other  federal  statutes  and  regulations  relating  to 
animals  and  experiments  involving  animals  and  adheres  to  princi¬ 
ples  stated  in  the  Guide  for  the  Care  and  Use  of  Laboratory  Animals 
(NRC,  201 1 )  in  facilities  that  are  fully  accredited  by  the  Association 
for  the  Assessment  and  Accreditation  of  Laboratory  Animal  Care, 
International. 

Citations  of  commercial  organizations  or  trade  names  in  this 
report  do  not  constitute  an  official  Department  of  the  Army 
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Appendix  A.  Supplementary  data 

Supplementary  data  associated  with  this  article  can  be  found,  in 
the  online  version,  at  http://dx.doi.Org/10.1016/j.bbr.2015.05.038 
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